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Abstract 

We define and calculate the discharge mode for a Schwarzschild black hole in massive 
electrodynamics. For small photon mass, the discharge mode describes the decay of the 
electric field of a charged star collapsing into a black hole. We argue that a similar 
"discharge of mass" occurs in massive gravity and leads to a strange process of black 
hole disappearance. 



1 Introduction and Estimates 

In massive electrodynamics, a charged star collapsing into a black hole must lose its electric 
field - the black hole must discharge. For a small photon mass m, the rate of discharge can 
be estimated as follows. One treats the mass term m^0, where cf) q/r is the electrostatic 
potential, as the density of the screening charge ps = —rri^cj)/ (47r). One further assumes that 
the screening charge moves into the black hole at the speed of light. Then q ^ Attt^Ps = 
—qm'^rg, where Vg is the Schwarzschild radius. This gives an exponential discharge, q oc e~'^^, 
with the decay rate 7 ~ m^r^,. In §2 we show that this estimate correctly describes the 
(intermediate asymptotic) decay of the electric field near the black hole in the limit of small 
photon mass. 

Now consider a black hole in massive gravity. The gravitational field of the black hole 

is screened by a negative energy density ps ~ —m^Adj,, where m is the graviton mass and 
Mp is the Planck mass. We argue in §3 that this negative energy must be accreted onto 
the black hole. Assuming, like we did in massive electrodynamics, that the screening energy 
accretes onto the black hole at the speed of light, we get a decreasing black hole mass: 
M ~ —r'^gW^Mp. Thus, the black hole loses mass and gradually disappears. 

The black hole disappearance is a weird prediction of massive gravity, as it seems to 
make possible the following scenario: (i) there was a star in an asymptotically Minkowski 
space-time, (ii) the star collapses into a black hole, (iii) the black hole disappeares leaving 
behind just the Minkowski space-time. 

To be clear, our results are inconsequential for real astrophysical black holes. If we live in 
an asymptotically flat universe with a massive graviton, the graviton mass must be smaller 
than the Hubble constant and the disappearance time much longer than the Hubble time. 
In fact, the phenomenological continuity of m — > limit is a quite satisfactory outcome. But 



the very possibility of disappearance (from the asymptotically Minkowski, eternal universe) 
seems to be a badness of the massive gravity. 

What makes the above scenario impossible in General Relativity is the conservation of 
the asymptotically defined ADM mass. In massive gravity the ADM mass vanishes - this 
is the characteristic property of the theory. And while an alternative globally conserved 
quantity can be defined in massive gravity (§3.5), global charges are not conserved in the 
presence of black holes. ^ 



2 Black hole discharge in massive electrodynamics 

In this section, we first formulate the Einstein-Proca theory and argue that charged black 
holes must discharge (§2.1). We next study the decay of the electric field and define the 
discharge mode (§2.2). Finally, we find quasi- stationary time-dependent solutions in the 
limit of small photon mass and identify the discharge mode (§2.3). Our massive gravity 
calculation will closely follow the calculation of §2.3. 



2.1 Einstein-Proca theory. Singularity of static charged black hole 
solutions 

The field of charged stars and black holes in massive electrodynamics is governed by the 
Proca equation 



d.i^g^'^g^'^F^p) + m^y^^A^ = A^i^J^', (1) 
where F^^ = d^j^Aj^ — d^A^, and is the electric current. The Einstein equation is 

n =T 4- T^™) (2) 

where T^^, is the sress-energy tensor of matter and T^IT'' is the stress-energy tensor of the 
Proca field 

= T^'J + m\A,A, - ]^g,,g-f'A^Ap), (3) 

where t'^^J is the stress tensor of the Maxwell theory. For black holes and outside stars 
= 0, TJ^ = 0. Here, we only consider the spherically symmetric case where the metric 
is diagonal in {t, r, 6, (p) coordinates. The spherical symmetry also implies that the only 
non- vanishing component of Ffj_,y is Ftr- 



^After submitting the first version of this paper, we were informed of the work [1] on black holes in 
bimetric gravity. There, it has been argued that asymptotically flat static black holes are non-singular only 
when the two metrics coincide with each other, and hence, with the Schwarzschild solution of the Einstein 
gravity. Thinking of massive gravity as a bimetric theory in the limit where the Newton's constant G of the 
second metric vanishes, the gravitational radius with respect to this metric is always zero, leading to the full 
disappearance of black holes in massive gravity. For any finite G, we expect the "mass discharge" of black 
holes in bimetric gravity to continue until the gravitational radii with respect to the two metrics coincide, 
at which point the time-dependent black hole solutions settle to the regular static solutions. 
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For static solutions, the only non-zero component of the Proca field is A^. The Ar 
component vanishes, because the ii — r component of equation (1) reads 

dtiV^gg'^g'^'F,,) + m^v^/'-A = 47^y=^J^ (4) 

which forces A^ — Q when dt — Q and J'' = 0. The ji — t component of (1) then becomes 

driV^gg^'g^'drA,) + m^^gg^'A, = Att^J^ (5) 

and gives A^. Solutions vanishing at spatial infinity have the expected Yukawa behavior at 
large r: 

At Ki when r » Vg. (6) 

r 

Defining the electric charge via the Gauss's law, this shows that localized sources are screened 
by the opposite charge density —m'^A^/ATr, carried by the Proca field. 

This is how one can calculate the electric field of a charged star. But static charged black 
holes do not exist, as they would be singular at horizon. To demonstrate this, consider the 
invariant product g'^'^A^^A,^, which is an observable in Proca theory. The product diverges 
at horizon because it reduces to g^^A'f and At ^ ai horizon (if At were to vanish both at 
horizon and at infinity, eq.(5) would give At{r) = 0, corresponding to an uncharged black 
hole). 

The horizon singularity of the static solution indicates that newly formed charged black 
holes would get rid of their electric field hair. The corresponding time-dependent solutions 
will be studied in the next section. 



2.2 Non-singular time-dependent black holes and the discharge mode 

For simplicity, we assume that the charge is small. Then we can use the Schwarzschild metric 

ds^ = (1 - '^)df - (1 - "^r'dr^ - rHn\ (7) 

The vacuum Proca equation for spherically symmetrical field can be transformed to an 
equation for the physical electric field E = —F'^ 12 (the field is, of course, radial) 

dlE - (1 - ^)5. ((1 - 'j:)]rM^'E)) + m\l - 'f)E = 0, (8) 

which can be brought to another useful form: 

{d^-dl + V)^^0, (9) 

where ^ — rE, p — r + Vg ln{r/rg — 1) is the tortoise coordinate, and the effective potential 
is 



3 



One can now set up an initial value problem by specifying the initial \E' and \E' and 
integrate eq.(9) forward in time. Since V is non-singular for all — oo < p < oo, an initially 
non-singular \l/ will remain so at all later times. Numerical experiments with various initial 
conditions confirm the anticipated decay of any initial electric field. Moreover, for small 
photon mass the decaying solutions asymptote to a mode that is best described by the 
intuitive picture of a quasi-stationary accretion of the screening charges onto the black hole: 
the discharge mode. 

We formally define the discharge mode by postulating the exponential in time decay 
of the field \E'(t,p) = e~''^'^{p). The field is in-going at horizon and decreasing at spatial 
infinity: 

i-i^/ + l' + V)^ = 0, (11) 
_ = -^^, p^-oo, (12) 

^-^0, p^+oo. (13) 

This eigenvalue problem is solved numerically in Appendix A, where we also show that our 
previous estimate of the discharge rate happens to be asymptotically exact: 7 = m^r^ in the 
limit mrg — 0. 

We emphasize that the discharge mode is not a true long-time asymptote. For t ^ 
(m^r^)"^, the electric field will become oscillatory with an algebraic decay t~^^^ sinmt, as 
shown in [2]. Nor does the discharge mode uniquely describe the time evolution of the 
electric field outside a charged star which collapses into a black hole. The evolution will 
depend on how exactly the collapse occurs, and the electric field will be given by a linear 
superposition of the decaying initial field and the field emitted by the charges moving into 
the black hole. If the charges move into the black hole fast enough, it seems reasonable to 
assume that the radiation will be mostly beamed into the black hole. The outside field is 
then probably dominated by the discharge mode. But in any case, the decay of the electric 
field cannot occur faster than the decay of the discharge mode. So to be precise, for small 
m, the intermediate-asymptotic decay of the near-hole electric field occurs in the discharge 
mode. 



2.3 Discharge mode as a quasi-stationary solution 

In the limit of small mrg, the quasi-stationary character of the discharge mode becomes 
sharper. In this section we show that there exists a two-parameter family of quasi-stationary 
solutions characterized by the initial charge of the black hole q and the discharge rate q. The 
previously defined discharge mode corresponds to the unique q/q ratio for which the horizon is 
non-singular. The approximate method used in this section helps us extract analogous quasi- 
stationary solutions in massive gravity, where on the one hand, time-dependent simulations 
seem impossible, and on the other, mrg ^ 1 is a reasonable assumption. 

Let us return to the vacuum Einstein-Proca equations (1,2). The first observation to 
make is that when mvg <^ 1, the expected discharge rate q ~ m^rgq is so small that the 
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induced time-dependence of the electric field and the metric are negligibly small as compared 
to the r-dependence. The only effect of time-dependence is to allow Aj. ^ 0. This can be 
seen from the expression for the screening charge current JJ = —m^J^j^.'K^ which gives 

^ = m^^/M,|^^^ . (14) 



The second observation is that in the near-hole region we have r ^ m^^, and as long 
as the fields are regular, the terms proportional to m in eqs.(l) and (2) are negligible. We 
therefore recover the time-independent Einstein- Maxwell theory, but now with an additional 
'gauge' condition on 

d,{^gg^^A,) = 0, (15) 

which is obtained by taking the divergence of eq.(l). This equation happens to be indepen- 
dent of m, and in our approximation it becomes an equation for Ar 

a.(v/=^/M,) = 0. (16) 

This is just the condition of stationary flow on J^. 

Assuming the Schwarzschild metric (7), we solve (16) 

A = (17) 

r[r — Tg) 

where c is a new integration constant related through eq.(14) to the discharge rate c = 
—q/w?. As for At, we can use the Einstein-Maxwell result At — q/r. Thus, we have a 
two-parameter family of quasi-stationary solutions characterized by the charge q and the 
discharge rate q. 

Requiring the horizon to be regular uniquely fixes the discharge rate. For instance, the 
regularity of the norm squared g^^A^^A^ = g^^A^ + g^^A^ gives 

q = ±m'^rgq, (18) 

which agrees with the previous calculation of the discharge mode in the limit mrg 0. 
The sign is undetermined by our equations, which are invariant under the time-reversal, but 
physical considerations select the negative sign.^ 



3 Black hole disappearance in massive gravity 

In this section, we first introduce Fierz-Pauli massive gravity, its nonlinear completions, and 
some of their relevant properties (§3.1). We next consider asymptotically flat spacetimes 
and argue that black holes must disappear (§3.2). In §3.3, we find quasi-stationary solutions, 
and show that unlike massive electrodynamics there is no unique disappearance rate. We 
next examine the implications of our findings in the actual problem of star collapse (§3.4). 
Finally in §3.5, we discuss the concept of global energy in massive gravity, and its failure in 
the presence of black holes. 

^The above calculation can easily be generalized to the case where the back-reaction of charge on geometry 
is not negligible. The near-hole metric will then be the Reissner-Nordstrom solution and Vg in (18) will be 
replaced by r+, the radius of the outer horizon. 
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3.1 Fierz-Pauli massive gravity 



By Fierz-Pauli massive gravity (FP), we denote a class of theories described by the Einstein- 
Hilbert action plus the (non-linear version of) Fierz-Pauli action: 

S = Seh + Sfp, (19) 
Seh = -U d'xy^R, (20) 
Spp = m^ J d^x^U. (21) 

Here we have set SnG = 1, and f/ is a Lorentz-invariant potential defined using a flat 
reference metric rjab- For small hah = gab — Vab it reduces to the Fierz-Pauli mass term [3] 

U^'^ = l{h'-hl). (22) 

where h = rj'^^hab- The Lorentz-invariance of U can be enforced by requiring that f/ be a 
symmetric function of the eigenvalues of the matrix = g°''^rii,c. 

The Fierz-Pauli term (22) is special because it is the unique quadratic expression that 
is linear in the perturbations of the ADM lapse function 6N = N — 1. So as in general 
relativity, remains a Lagrange multiplier and the theory describes five dynamical degrees 
of freedom at quadratic level. Rather obviously in the hindsight, this will no longer be the 
case beyond the quadratic level for a generic U, as Boulware and Deser pointed out [4]. This 
results in a sixth degree of freedom which they showed to be a ghost. This problem has 
been solved by de Rham, Gabadadze, and ToUey [5] who found a particular two-parameter 
family of potentials which propagates just five degrees of freedom at all orders. A concise 
representation of the family is [6] 

[/ = ^ \a\b + C2 A„AbAc + C3A0A1A2A3, (23) 

where the sums are over all all-distinct pairs and triples of indices, and A^ are the four 
eigenvalues of the matrix 

5» - ^/^. (24) 

It is easy to understand why this two-parameter family, which we henceforth refer to as 
FP2, is special. Restricting to metrics with zero shift vector, the lapse appears only in 
Ao = 1 — A^^^. Apparently, the expression (23) is the only symmetric combination of Aa, 
such that y/—gUFP2 has zero cosmological constant, and is linear in A^. The full proof of the 
absence of Boulware-Deser ghost is given by Hassan and Rosen [7] (see also [8] for a different 
approach).'^ 

As formulated above, the FP theory is not generally covariant. However, it can be made 
so by introducing 4 non-canonical scalar fields [14, 15, 16] to write as 

= Vbcg'^'d^rd,^". (25) 

■^Evidently, having the right number of degrees of freedom does not protect FP2 against other pathologies 
such as strong coupling, instability, or superluminality of fluctuations when the background metric deviates 
from Minkowski (see, e.g., [9, 10, 11, 12]). The asymptotically flat solutions considered here are not expected 
to be exceptional [13]. 
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In the so-called unitary gauge, one uses these scalar fields as coordinates, X" = 0"^, and 
recovers the original formulation of the FP theory. In the unitary gauge, the theory is still 
invariant under simultaneous reparametrization of Qat and rjab- For our purposes - studying 
spherically symmetrical stars and black holes - we transform to spherical coordinates. Then 

r]abdX''dX^ = dT^ - dR^ - R^dn^ . (26) 

If the space-time is spherically symmetric, the most general metric in this coordinate system 
is parametrized by four functions of T and R: 

ds^ = CdT^ - 2DdTdR - AdR^ - BR^d^f . (27) 

A comment on our notation: in the following we will frequently switch to the time 
variable t as measured by asymptotic observers, and the circumference defined radius r (the 
Schwarzschild variables). However, the numerical indices (0,1,2,3) are exclusively used to 
denote unitary-gauge variables in the spherical coordinates in the order (T, i?, 9, (p). 

3.2 Asymptotically flat spacetimes and the inevitable singularity of 
static black holes 

As in massive electrodynamics, it is natural to first look for static star and black hole solutions 
in massive gravity. The procedure is briefiy outlined in appendix B, and it was fully pursued 
in [17] with the conclusion that there exist acceptable star solutions but black holes are 
generically singular at horizon. However, again as in the case of massive electrodynamics, it 
is easy to exhibit the horizon singularity of static black holes without knowing the explicit 
solution. The horizon singularity now indicates that a star collapsing into a black hole must 
get rid of its very gravitational field and hence disappear.^ 

The horizon singularity can be shown as follows. Consider the most general static spher- 
ically symmetric vacuum solution of massive gravity. The unitary-gauge metric is given by 
(27). Requiring the space-time to be asymptotically flat then forces the metric to be diagonal 
[17], i.e. goi = D = 0: 

First note that the Ricci tensor of a time-independent metric of the form (27) 
satisfies the identity [18} 

5'01-Roo - fl'oo-Roi = 0. (28) 

This imposes the following purely algebraic constraint on the unitary-gauge metric 
components via the vacuum massive gravity equations: 

QqiTqq - gooToi = , (29) 

where Tab is the stress-energy tensor of Spp- 

^Sergei Dubovsky pointed out that the horizon singularity of our static black hole solution [17] might 
signal accrertion. like it does for fluids. If one tries to find a static fluid surrounding a black hole, one gets 
a singularity at horizon; in reality fluids are accreted. 
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On the other hand, for any potential U that is a symmetric function of the eigen- 
values of the matrix = g"''^T]cb, one can verify that Tqi = kqoi, where the 
proportionality coefficient n is a non- singular function at g^i = (see footnote 
13). It follows that eq. (29) divides the solutions into two branches 

(i) 9oi = 0, 

(ii) Too = i^goo- 

Now consider large radii R ^ oo, where the metric perturbations hab = gab — Vab 
are small by the asymptotic flatness assumption. For any non-linear completion 
of the Fierz-Pauli mass term /i^j, — h"^ , we then have 

Tm = \ni\oi{l + O{h)), (30) 

giving k = m? /2 at large R. This excludes the branch (ii) as it leads to a fi- 
nite asymptotic value forT^Q. This branch is where the Schwarzschild-de Sitter 
solutions of fl8, 19] are realized (see also f20j in that context). 

On the other hand, whenever the unitary-gauge metric is diagonal, horizons will be 
physically singular in massive gravity [21]. This is because the inverse unitary-gauge metric 
components are scalar quantities in massive gravity: g""^ = g^'^dn(j)"'d^(f)''. When g"'^ is diago- 
nal, its {00} component will be singular at horizon and this singularity will be reflected in 
the action and the stress-energy tensor via Aq = 1 — \^g^. To have finite g*^^ at horizon, 
one necessarily needs gQi ^ which, by the above arguments, leads to time-dependence and 
non-zero energy fiux Tq. ^ 

In view of black hole discharge in massive electrodynamics, the time-dependence of black 
holes in massive gravity is in fact naturally expected. The linearized field of a point source is 
known (and shown in §3.3.1) to exhibit the Yukawa decay hab oc exp(— mr) in massive gravity. 
It follows that the invariant ADM mass of any localized system is zero. Therefore, the mass 
term, taken to the right hand side of the Einstein equation, can be thought of as the stress 
tensor of screening matter with negative energy that surrounds and degravitates gravitational 
sources. When a black hole forms, the screening matter fiows inside and diminishes the black 
hole mass.^ 

The rate of this process is estimated by requiring the near horizon energy density Tq ~ 
— accrete with the speed of light through an area of order r^: 

f-g mVJ, (32) 

^Another way to see why goi 7^ is a necessary condition for regularity of the horizon is to realize that 
when goi = 0, the unitary-gauge time variable T ~ (jp coincides with the proper time of the asymptotic 
observers t. But t is a singular variable at horizon and therefore cjp — t i& singular [22]. This inevitable 
singularity relies just on the existence of a horizon, but if the Schwarzschild geometry is indeed recovered at 
short distances a la Vainshtein, we know how exactly t diverges at horizon: in terms of the advanced time v 
and circumference defined radius r which are regular variables 

t — V — p — V — [r + Tg ln(r/rg — 1)]. (31) 
^The classical instability of black holes in massive gravity has been conjectured before by Gia Dvali [23]. 
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where for a mass M black hole Vg = M/Att in our unites. Note that when mrg <^ 1 the 
associated scale r — l/m?rg is much longer than other length scales in the problem. We will 
use this fact to find approximate time-dependent solutions in the next section. 

Assuming that the time-dependence is mild, of order (32), we can estimate the resulting 
value of the off-diagonal metric component g^i. Note first that the identity (28) is violated 
in the time-dependent case by terms of order 

QqiRoo - 9mRo\ ~ ^o^iS'oo ~ Tg/r'^ ~ m^ryp, (33) 

where we used the estimate (32) in the last step. Next replace Rab by Tab via the Einstein 
equation, as we did to find the constraint (29), and use the same arguments to find 

501(^00 - f^Qoo) ~ mVJ/r^ (34) 

Since (701 = as r — >■ cxd and k ~ m^, we get qqi ~ r|/r^. Thus, we expect the diagonal 
static solution to be a good approximate solution everywhere except very close to the horizon 
where Qqi becomes of order unity and makes the inverse metric g""^ finite. This expectation 
will be confirmed by the calculation of the next section. 

3.3 Non-singular time-dependent black holes 

It is relatively easy to find exact time-dependent solutions in the linear regime, but this is 
not the case at non-linear level. We therefore break the problem into three parts: first, we 
find the exact linear solutions and show that they are characterized by two parameters, the 
mass Tg and the disappearance rate fg. Then we use the mr-g <^ 1 approximation to find 
quasi-stationary near-hole solutions, again parametrized by the mass and rate. Finally, we 
match the linear and the near-hole solutions. We will see that unlike massive electrodynamics 
where all but a unique ratio q/q lead to a singularity at the horizon, here all quasi- stationary 
solutions are regular. The analogous case of the accretion of superluminal fluids and possible 
explanations of this result are discussed in the following section §3.4. 

3.3.1 Linear solutions 

To study the linearized system, we parametrize the metric as 

ds'^ = (1 + c)dT'^ - 2d dTdR - (1 + a)dR^ - (1 + b)R^dn'^, (35) 

where a,b,c,d are infinitesimal functions of T and R. The linearized Einstein equation in 
Cartesian coordinates reads 

Dhab - dadchl - dbdjil + dadbh - rjabOh + rjabdcddh"'^ = m'^{r]abh - hab), (36) 

and its divergence yields the condition dah'^ = dbh. Since the equations are first order in 
perturbations hab, we can replace in them (T, R) with (t, r) which coincide at zeroth order. 
Moreover, the solutions can be expanded in the exponential basis 

a, b,c,d(x e-^^^K (37) 
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The system of equations can then be brought into the form 

A2 



a'r + 2a - 26 = —^^m^r'^ (a - 2X%) , (38) 
1 — 2A^ 

h'r + h-a= ^^^_^^^,^ mV(a - 2\%), (39) 



where ' = d/dr, and c and d are given in terms of a, 6 by 

cV = 2(1 + A^)(a-6-6V), (40) 
2A 

= — ia-h-h'r). (41) 
mr 

Thus, for any fixed rate |A| < l/-\/2 there exists a one-parameter family of decaying as a 
function of r solutions. The parameter characterizes the mass of the gravitating body. 

Note that in agreement with the arguments of the previous section, g^i vanishes in the 
static limit A — > 0. For an accretion rate of order (32), we have A = amrg ^ 1 (where 
we introduced a as a so far undetermined order-unity parameter characterizing the rate Vg). 
The solution can therefore be approximated by the A = vDVZ solution [24]: 

c = --^e-'"^ a = --^e-™'-(l + mr), 6 = ^^^e-'""(l + mr + mV), (42) 
3r 3m^r'^ Sm'^r'^ 

except that Vg adiabatically changes with time as rg{t) ~ rg(0) exp(— am^r^t), and 

4ar^ 

d = —^e-'^'il + mr). (43) 

We see that, apart from allowing small non-zero qqi, the time-dependence is inconsequential 
in the long-distance linear field.'' In FP2, this linear solution is valid for r S> ry = {rg/m?Y^^ . 



3.3.2 Near-hole solutions 

Similar to massive electrodynamics, we can develop an approximate method to find the short 
distance (r ^ ry) quasi-stationary solutions in the limit mr^ ^ 1. This approximation is in 
fact a generalization of the Vainshtein's original idea [27], which we now review and extend 
to our case. 

Aiming for a static solution, Vainshtein started from a diagonal metric ansatz 

ds^ = CdT'^ - AdR^ - BR^dn"^, (44) 

^We also see that in the context of massive gravity the no-hair theorem of [25] and in general the 
singularity of static black holes cannot be interpreted as a dramatic difference between the field of stars and 
black holes. This can be understood from the fact that the linear static field of stars and black holes is 
uniquely determined by the mass r^. To lose its hair, the black hole must actually lose its entire mass, which 
is done only very slowly. In particular, we think that contrary to the claims of [26], the response of stars 
and black holes of the same mass to the field of distant objects is almost indistinguishable. 
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where A, B, C are functions of R, and are determined using the Einstein equations 



G^, = Ti^. (45) 

He noticed that is proportional to m? and therefore can be ignored at short radii, as long 
as the metric coefficients remain finite. Therefore, the mass-less Einstein theory of gravity 
is recovered in this limit, implying that the metric (44) must be a reparametrization of the 
Schwarzschild metric 

ds'' = (1 - '^)de - (1 - ^)-Mr2 - rHn\ (46) 

The reparametrization is determined from the covariant divergence of eq.(45): 

V,T!^ = 0, (47) 

which can be thought of as a gauge condition on the metric (44). Note that m drops out of 
this equation, so it is a non-trivial constraint even in the zero-graviton-mass limit. 

Requirement of asymptotic fiatness {gab = Vab as i? — )■ oo) then fixes T = t, and the only 
non-trivial component of eq.(47) (the u = 1 component) serves as an equation for R{r), in 
terms of which A,B,C are given by 

C = l-^, A = 1/(CR'^), B = ryR^. (48) 
r 

Vainshtein realized that while there is no such reparametrization of the Schwarzschild metric 
at the linearized level (the well-known vDVZ discontinuity), linearization becomes inadequate 
at r ~ ry. He showed that at the non-linear level there exists a finite solution valid for 
<^ r ^ ry, which a posteriori justifies neglecting T^^ from the Einstein equation. One still 
needs to check whether this solution matches the vDVZ solution for r ^ ry, where the mass 
term cannot be ignored anymore. This is a non-trivial check, but there exists a sub-family of 
FP2 for which the answer is positive [17, 28] (see [29] for earlier works, and [30] for a related 
discussion). 

Returning to the time-dependent problem, we expect that, as in massive electrodynamics 
where for small mrg the only relevant effect of the time-dependence was to excite Ar, here 
the time-dependence excites qqi. Once this is taken into account by using the metric ansatz 

ds^ = CdT^ - 2DdTdR - AdR^ - BR^dn^, (49) 

we can ignore time derivatives (quasi-stationary approximation) and repeat Vainshtein's 
procedure. Neglecting the stress tensor of the FP action from the Einstein equation implies 
that the metric is a special reparametrization of the Schwarzschild metric which satisfies the 
gauge condition V ^Tj^ = 0. 

In the quasi-stationary approximation, the metric (49) can be diagonalized and expressed 
in terms of t and r by the coordinate transformations 

R = B-^/'^r, (50) 

T = t+ 1 ^^dr, (51) 
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where 7 = dR/dr = B ^/^(l — B'r/2B). In this coordinate system the metric looks like 

ds^ = Cdf -i\A + D^/C)dr^ - r^dn\ (52) 

Requiring this to match the Schwarzschild metric (46) gives 

C = l-rg/r, A = C-\-f~^ -D^). (53) 

Hence, there are two unknowns D and B, characterizing the gauge transformation from 
the unitary frame to the Schwarzschild frame. The task is to solve for them using the 
u = 0,1 components of V^Tj^ = on a fixed Schwarzschild geometry, and subject to ap- 
propriate boundary conditions. The resulting quasi-stationary solutions are expected to be 
parametrized by the mass and the disappearance rate Vg = Q;m^r| (which will be the 
integration constant of the i/ = equation). One must then ask for what values of a the 
solution is regular at horizon.^ 

Note that the above procedure is equivalent to solving the equations of motion for the 
scalar fields 0"^ on a fixed Schwarzschild background (as one usually treats accretion prob- 
lems [31]). By spherical symmetry 0* = Rn^ = B^^^'^x\ and the stationary approximation 
corresponds to looking for solutions of the form 0° = t + (f{r) [which is related to B and 
D by eq.(51)]. This is why there are only two non-trivial equations; they are related to the 
stress-energy conservation via the Bianchi identity 

= -^d,r^^, (54) 



-9 

which is particularly simple in the unitary gauge where (9b0" = 6^. Neglecting the back- 
reaction on geometry is justified for small graviton mass if we find a regular solution. 

The details of the calculation for FP2 is given in appendix C. In summary, we obtain an 
algebraic equation for A2 = 1 — -B"^/^: 

3r / 

^(1 - ^) + (C2 + C3A2)(1 + (1 - X2)C) =[C+ ^ j (1 + C2X2 + 2(C2 + C3A2)), (55) 

where 02,3 are related to the parameters of FP2 action (23), and P = 1 + 2C2A2 + C3AI is a 
positive function when < C3, the range where stable numerical time-independent solutions 
exist. Given a and Vg, eq.(55) is solved for A2, using which D is given by ^ 



2^A\ -1/2 



D = ^(r'^C)[C+^^) . (56) 

^In the language of footnote 5, one seeks solutions on which 0" interpolates between (fP — t a,t r — 00 
and = t + p at r = [22]. 

^The uniqueness of the solution once and fg are fixed is a sign of the absence of Boulware-Deser 
ghost in FP2. Among the five degrees of freedom of a massive graviton only the scalar mode is dynamical 
in the spherically symmetric problem, and can participate in the accretion (analogous to the longitudinal 
polarization of the massive photon and the sound mode of fiuids) . After the large-r asymptotic condition 
is fixed by the mass Vg (charge q for photon, asymptotic density for fiuids), only one additional integration 
constant is needed to fully parametrize the stationary accretion. In generic FP, on the other hand, the 
analogue of equation (55) for B is of the second order, which requires two extra integration constants. These 
correspond to the emergence of a second scalar degree of freedom, which is the Boulware-Deser ghost. 
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Now we can estimate how much the time-dependent solutions deviate from the static 
solution. Assuming a ~ 1, equation (56) gives a small off-diagonal metric component at 
r > rg-. 



and eq.(55) reduces to its static a = limit, giving A2 = l/^/c^. Therefore, as anticipated 
in §3.2 the effects of time-dependence are felt merely within a distant of a few Vg from the 
horizon, beyond which the static solution of [17] remains accurate. In particular, matching 
to the linearized solution (studied in appendix D) does not impose any constraint. This 
should not be very surprising in view of the fluid accretion problem. The critical radius that 
determines the accretion rate for fluids is close to the horizon for relativistic fluids. This is 
also what one expects from the example of massive electrodynamics where regularity of the 
horizon itself determines the accretion rate. 

What is the actual rate of accretion r^? In massive electrodynamics, the discharge 
solution has been selected from the two-parametric (g, q) family of quasi- stationary solutions 
by the horizon non-singularity. If massive gravity was fully analogous to electrodynamics, 
the rate Vg would also be uniquely determined by the horizon non-singularity. But this does 
not actually happen. For r — )■ (and C — )■ 0), eq.(56) gives 



and upon substitution in (51) gives 0° = T = t + ep + 0{a^^), which is regular at black 
hole horizon for any positive a (see footnote 8)}^ So in contrast to massive electrodynamics, 
any non-zero rate seems to result in a regular solution in massive gravity. This behavior, 
though peculiar, is not unprecedented. It is known [32] (and reviewed in appendix F) that 
the accretion rate of fluids with superluminal speed of sound onto black holes is not unique. 
In the next section, we examine the implications of this finding for the actual process of star 
collapse, where we argue that the rate is determined by the history of collapse. 

3.4 The fate of collapsing stars 

In classical physics, an outside observer never sees a fully formed black hole - collapsing 
stars are stuck forever at their gravitational radius with redshift increasing at their surface. 
In this section, we study black hole accretion from this perspective which, besides being 
more realistic, provides a useful thought laboratory to examine our analytic results. We first 
consider the accretion of fluids onto collapsing stars. We then give a related interpretation of 
black hole discharge in massive electrodynamics, and finally, turn to the problem of collapsing 
stars in massive gravity. 

^''if one only requires the regularity of g"'', then negative values of a are also admissible. In that case, the 
black hole excretes negative energy 'aether' and grows. 




(57) 





(58) 
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Figure 1: The density of static fluids (represented by the depth) always diverges at horizon. 
However, the total amount of fluid mass that can be stored at rest between the star surface and the 
horizon is infinite for < 1, while finite for > 1. 



3.4.1 Fluids 

Ideal fluids are characterized by their density, pressure, and velocity fields, respectively e,p, 
and u'^, and a certain equation of state that relates e and p. The fluid stress-energy tensor 
is 

%i, = {e + p)u^Uu - pg^,u- (59) 

Consider a spherically symmetric static solution on the Schwarzschild metric (we ignore the 
fluid back-reaction on geometry). The u = 1 component of the stress-energy conservation, 
"^i^V = 0' gives 

Assuming a simple p = c^e equation of state, we get the density field 

e = Boil - ^)-(^+- (61) 
r 

where Eq is the density at spatial infinity. We see that for subluminal and superluminal fluids 
alike the density and pressure diverge at r = r^, reflecting the fact that an infinite force is 
needed to hold a test particle at rest above the horizon. 

However the two cases, subluminal and superluminal, are known to be very different 
when one consideres the black hole accretion. The rate of accretion is uniquely determined 
by the critical Bondi solution in the subluminal case, < 1, because accreting solutions 
with slower rates are singular at horizon, while solutions with faster rates do not exist. On 
the other hand, in the superluminal case, when > 1, the critical accretion rate does not 
exist and all accreting solutions are regular at horizon (see appendix F). 

This difference is ultimately due to the different amounts of static fluid mass which can 
be stored near the black hole horizon. On Schwarzschild metric, there is a well-defined 
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notion of fluid mass since the i/ = component of the stress-energy conservation becomes 
9iJ-{V~~9%^) — imphes the conservation of 



M 



= J d'x^X- (62) 



At spatial infinity where the fluid is at rest and the metric is flat this reduces to J d^re. 
Hence, M is the fluid energy as measured by asymptotic observers. Now consider a high 
redshift star of radius Tq as a model for the black hole. Using the static solution (61), it is 
easy to see that when < 1 the fluid mass near the star surface diverges as ro — )■ r^, while 
it remains finite for > 1 (figure 1). 

Imagine next initiating the fluid accretion by contracting the star (classically, ro will 
pass the horizon only at t = cxo, no matter how fast the contraction takes place, so the 
accretion problem can be fully addressed in this framework). When < 1, the above- 
mentioned divergence of near-horizon fluid mass makes it possible to sustain an inflnitesimal 
rate of steady accretion by adiabatically contracting the star. The pressure and density in 
this regime are then approximately given by the singular static solution. This explains the 
singularity of those sub-critical accreting solutions. In contrast, because of the finiteness 
of near-horizon M in the > 1 case, it is impossible to sustain any steady accretion by 
adiabatic contraction of the star. Consequently, the accreting solutions of superluminal fluids 
can never be approximated by the singular static solution near the horizon. 

Without an explicit calculation of fluid accretion we cannot proceed much further. So in 
the following, besides filling in some of the details of the above picture, we use the results of 
appendix F to provide a more comprehensive description. 

Let us define x — r — rg in terms of which the t — r part of the near-horizon Schwarzschild 
geometry is given by the Rindler metric: 

ds^ — r~^xdt^ — TgX'^dx'^. (63) 

For a radially moving time- like curve near the horizon, we have 

^2 ^ rg^x\l - Tg^xe-^), (64) 

where the overdot denotes f-derivative, and e = Uq is the energy per unit rest mass as 
measured from infinity, e is constant for free-falling observers, but it vanishes as e = (x/r^)^/^ 
for static ones. Thus, the maximal value of near the horizon is |i;|max = ^/f^g- 

As a function of the position of the star surface xq = ro — Vg, the mass of static fiuid with 
density (61) that resides near the surface scales as 

xf^^^". (65) 

When < 1, this diverges in the limit xq ^ and, as explained above, infinitesimal rates 
of steady accretion can be obtained by adiabatically decreasing xq- 

M = -AnVgEXo, (66) 
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with e given by (61). As |xo| is increased to the point that this rate exceeds the critical 
Bondi rate, the fluid dynamics decouples from the motion of the surface and settles to the 
regular Bondi solution. 

On the other hand, when > 1 substituting the maximal value |xo|max = ^a/^g in eq- 
(66) gives zero as Xq — )■ 0. Thus, similar adiabatic contractions do not correspond to any 
steady accretion for superluminal fluids. This result may seem rather counterintuitive as it 
suggests that steady accretion is impossible when > 1. This is of course not the case. 
Consider again the near horizon contribution to M for a possibly moving fluid: 

M{xq) r-^ Aixrl I dx[{e + pYgC^x''^ - p], (67) 

where the expression in the square brackets is Tq for a fluid with Uq = e. Substituting the 
inflnitely redshifted e = (x/r^)^/^ for a static fluid, we see that the density of superluminal 
fluids is not singular enough and the integral converges as Xo — )■ 0. However, if the same 
fluid is allowed to fall in, with e 7^ at x = 0, even with a flnite e and p an infinite amount 
of mass can be deposited near the horizon. The flux 

Tq = {e + p)rge^x-^x, (68) 

would then be non-zero. The explicit solution of appendix F shows that any non-zero accre- 
tion rate of superluminal fluids can indeed be obtained by a speciflc choice of e at horizon. 
So different free-fall (e = const.) motions of the star surface at horizon lead to different rates 
of accretion, implying that unlike subluminal fluids the accretion of a superluminal fluid 
never decouples from the star collapse. Moreover, since even for inflnitesimal accretion rates 
the fluid is free-falling at horizon, it is expected that the pressure remains flnite there. 

Given that at any flnite t all of the accreted matter is still above the horizon, one may 
wonder what happens if the contraction of the star is paused after a long period At of steady 
accretion. Clearly the static solution must be recovered since the star surface is at a flnite 
redshift. The flniteness of the near horizon static fluid mass when > 1 then implies that, 
for large enough At, the star must excrete most of the fluid that has been accumulated near 
the surface back to larger radii. 

Finally, when = 1 the total mass of static fluid adjacent to the surface diverges loga- 
rithmically with Xq. Therefore, in the adiabatic approximation (66) even inflnitesimal rates 
require very fast motion (|xo| Xq, although with e oc x^^^) which makes this approxima- 
tion unreliable. A more detailed analysis is needed to show the singularity of sub-critical 
accreting solutions in this case. 

3.4.2 Massive electrodynamics 

Now let us interpret the black hole discharge in massive electrodynamics by considering a 
charged shell of radius rg approaching its gravitational radius Vg. We flrst need to determine 
the relation between the asymptotic fleld At = qe~"^^' /r and the actual charge Q on the shell 
as a function of tq. This can be obtained from the discontinuity of the electric fleld at the 
surface in the following way (a more detailed analysis can be found in [33]): 
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Ignoring back-reaction on geometry, and substituting At = qx/ r in the source-free 
Proca equation, we get for r > Tq 



T 



(69) 



For small rnVg, the decaying solution at infinity % = e goes into 

X — i + rn^rgX In x 



(70) 



near the horizon. On the other hand, inside the shell r < Tq, where Qrr = and 
Qtt ~ xo/rg = const.; we have 

X" = rn\. (71) 

Imposing regularity condition at the center and matching to the outside solution 
we get for r < Tq 

q sinh mr 

" ^ (72) 



ro mr 

The charge on the shell is then given by the discontinuity of proper electric field 
at ro 

+ r 1 

- qil 

'0 



Q 



gg^'g^Frtl'L ~ qm\ 



ln{rg/xo) + -{rg/xoy/^ 



(73) 



where we have kept both the contributions from outside and inside although the 
former is sub-dominant in the — >■ limit. 



To have a fixed asymptotic field (characterized by q), it is seen from eq.(73) that the total 
charge on the shell Q must be increased indefinitely as the surface redshift is taken to infinity. 
This is due to the divergence of the total amount of screening charge both inside and outside 
of the shell, respectively as and Inxo- In the realistic problem, as a star of fixed charge 
Q contracts to smaller radii, its asymptotic field diminishes since the surrounding screening 
charges accrete into the now deeper potential well around and inside the star. 

As in the case of subluminal fiuids the divergence of total static charge in the limit Xo ^ 
allows us to reproduce small rates of accretion by slowly contracting the shell: 

dQ . 1 2 5/2 -3/2 . 

q = ~Q^^o ^ ^qm rg' Xq xq. (74) 

However, here the screening charges do penetrate inside the shell and the slow motion of 
the surface does not imply that of the charges. In fact, the logarithmic divergence of charge 
density outside the shell makes this case analogous to fiuids with Cg = 1, and a detailed 
analysis is needed to show why discharge rates of less than m^r^ lead to a singular horizon. 
As the shell collapses faster, the outside field eventually decouples from the motion of the 
shell and decays via the discharge mode. 
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3.4.3 Massive gravity 

For static black holes in massive gravity, we show in appendix B.l that while pressure usually 
diverges at the horizon, the energy density always remains finite in FP2, and so does the 
total amount of static energy that can be stored outside the horizon. Therefore, as in the 
case of superluminal fiuids, to have any non-zero rate of accretion the surrounding screening 
matter should freely fall at horizon and hence the pressure is expected to remain finite. This 
seems to be the reason for the non-uniqueness of the disappearance rate. The rate will then 
depend on the motion of the collapsing star according to one of the following scenarios: 

i) If the total amount of negative energy inside static stars diverges with their redshift 
(similar to the divergence of total screening charge in the interior of the charged shell in 
massive electrodynamics), then by slowly contracting the star to higher redshifts we obtain 
arbitrary rates of accretion. This impfies that the gravitational radius Vg (which is approx- 
imately the actual radius of the star) shrinks as the redshift increases, and the perceived 
mass of the star decreases. In this case if the process of collapse is paused at some point the 
result will be a static star of smaller gravitational radius. 

ii) Otherwise, as in the case of impenetrable stars surrounded by superluminal fiuids, 
non-zero accretion rates are obtained when the surface freely falls which allows for the accu- 
mulation of an arbitrary amount of free-falling screening matter. As before, the gravitational 
radius (and the real radius) of the system shrinks in this process. However, if we decide to 
pause the collapse, the excess of screening matter must return to larger radii, resulting in an 
excreting star with ever increasing surface redshift. 

We were unable to find high redshift static star solutions in FP2 to decide between the two 
scenarios. Nevertheless, in both CclSGS clS the star shrinks to very small radii the quantum 
mechanical effects become important and the above description in terms of stars breaks 
down. The resulting small mass black hole will probably evaporate via Hawking radiation. 

3.5 A substitute for the ADM mass 

Let us conclude by a discussion of global energy in asymptotically flat solutions of massive 
gravity. We have already seen that the invariant ADM mass of a localized system trivially 
vanishes in massive gravity because of the Yukawa screening. Equivalently, if one defines an 
ordinarily conserved pseudo-tensor of stress and energy, say TJ^, the conserved total energy- 
momentum four-vector of the system which is obtained by integrating T° over the whole 
space is always zero. Intuitively, the negative contribution of Spp cancels the energy and 
momentum of the 'matter content', by which we mean everything except Spp. Of course, on a 
curved space-time these different contributions are not separately conserved, and seemingly, 
nothing forbids the growth or disappearance of the matter content at the expense of the 
growth and disappearance of the FP content. However, there exist four Noether charges, 
associated with four global symmetries of the Fierz-PauH theory, which as we will see are 
closely related to the negative energy-momentum of Spp- As such, they are also a well- 
defined measure of the matter-content of the theory. The conservation of these charges, 
then, constrain processes that involve stars, but they can be violated in the presence of 
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black holes. 

The prescription (25) to make Spp covariant, also introduces four global symmetries 
under constant shifts of the scalar fields 0". Therefore, there are four conserved Noether 
currents obtained by varying Spp with respect to = 9^0" 

and associated to them, there will be four conserved charges Qa- Moreover, the form of the 
Lagrangian U {g^^ e'^^e^j^'qcb) allows us to relate to the stress-tensor of Spp 

J^^eliTI^ + m'Si^U), (76) 



a 
b-> 



where T^y = 2{—g)~^/'^6SFp/5g^'', and is the inverse of e^. In the unitary gauge = 5, 
and therefore, is invertible as long as this gauge exists. Consider now a localized material 
system whose center of mass is at rest; Qq has the following properties: 

i) When the matter distribution is so dilute that the perturbations of the unitary-gauge 
metric are infinitesimal (/i«5 < 1), we have = + 0{hlf,). It follows that Qo = -M, 
where M = J d^rT^ is the total energy of the matter distribution. This can be seen from the 
fact that when linearization is possible, 1^ is simply the total stress-energy tensor appearing 
on the right hand side of the Einstein equation, so we have 



M, 



ADM 



= J dh%'o = J dh{T^ + 7^) = 0. (77) 



The condition hab < 1 is satisfied as long as G7^ < and GM'^I^T^^^^ < 1. 

ii) Inside the Vainshtein region, the gravitational field of a star is the same as in the 
Einstein theory. For a star at finite redshift one can show that M = —Qo, where M — AnVg 
is the Schwarzschild mass of the star as determined from its Einsteinian field: 

Consider a sphere of radius ri where Vg <^ ri <^ vy = {rg/rn^y^^. The space-time 
is already nearly flat at ri, so the ADM mass inside the sphere is well-defined 
and equal to M. The ADM mass inside a much larger sphere of radius r ^ oo 
is zero, so we must have 

/•oo 

/ dhT^ = -M. (78) 

J n 

In appendix E we show that 



/ dh4= / d^rT^. (79) 



Moreover, the assumption of finite redshift ensures that = 0{m'^) and the 



contribution of the region r < ri to Qo is negligible, because J^rf <C J^r 
0{rg). This fact together with equations (78) and (79), implies M — —Qq. 



^^We thank Gregory Gabadadze for discussions and collaboration that lead to the identification of these 
charges. 
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iii) If the redshift of the star diverges as inverse powers of mrg, large amounts of screening 
matter (presumably with negative energy) can be accumulated in the region r ~ r^, in which 
case the Schwarzschild mass of the star decreases and deviates from Qq. 

iv) When a black hole forms and evaporates conservation of global charges such as Qo 
can be violated. Therefore, the vanishing of the asymptotically defined Madm in massive 
gravity allows the post black hole matter content to be less than the initial one. 
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A Discharge mode in massive electrodynamics 

We numerically solve the eigenvalue problem 

i-i^/ + ^' + V)'b = ^^ ^ = _^vi>|_^, * = 0|+oo. (A.l) 

We take ^(pi) = e~'^^^, ^(Pi) ~ —76"'"'''% where pi is an arbitrary (but negative and large 
in absolute value) number. We then integrate eq.(A.l) from pi to large positive p. We choose 
7 so as to get — )■ 0, p ^ +00. This gives the decay rates 7. 

The resulting decay rate, for mrg < 0.5, to better than 3% accuracy, is given by a fitting 
formula 

7 fti (A.2) 

The asymptotic 

7 — >■ m^rg, mVg — >■ (A. 3) 

can be shown to be exact. One first calculates the (singular at horizon) stationary electric 
field from the Proca equation (8) 



dr \ J- dr 
It follows that 

V + (1 --)-* = -"iV / drE. (A.5) 
dp r r Jj. 

For small m, the right-hand side of eq.(A.5) can be approximated by the zero-mass solution 
E = q/r^, giving 

^ + (l-^)-* = -mV (A.6) 
dp r r 
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For r close to the horizon, r — ^ r^, we get 

— = -m^q. (A.7) 
dp 

The decay mode, for r — ^ r^, is \1/ = Ce~'^'^ where C is some constant, so that 

- = -0'*. (A.8) 

On the other hand, for negative p in the interval 1 ^ \p\ ^ 7^^, the decay mode and the 
zero-mass static mode are approximately equal to the Coulomb field at horizon, \1/ ^ q/tq, 
and comparing eqs.(A.7, A.8) we get 7 = m^Tg. 



B Static spherically symmetric field 

In this appendix we outline the method used to find the static solution and a few relevant 
results of [17]. Starting from the unitary gauge, the asymptotically flat metric can be written 
as 

ds^ = e^'dt^ - e^dR^ - R^e^d^^. (B. 1) 

We then change the radial coordinate R to the circumference defined r 

ds^ = e'^dt^ - e^dr^ - r^dQ^. (B.2) 

The scalar fields in the new coordinates become 

/ = t, 0^=re-'^/V, (B.3) 

where is the unit radial vector. 

As independent vacuum equations we use two Einstein equations and the stress-energy 
conservation 

= r-2(l - e"^) + r-'e-^X' = T°, (B.4) 
Gl = r~\l - e-^) - r-^e-\' = (B.5) 

r^' = iz.'(T°-T/) + ^(r|-r,^), (b.6) 

where prime denotes the r-derivative. Inside stars one adds the matter stress-energy tensor 
to the right hand side of eqs.(B.4,B.5), and separately imposes its conservation JT^ = 0. 

The stress-energy tensor of FP2 is derived from eq.(23): 

T° = -m\Xi + 2X2 + C2(2AiA2 + A^) + cgAiA^), (B.7) 

Tl = -m^{Xo + 2X2 + C2(2AoA2 + A^) + csAqA^), (B.8) 
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T| = -m^{Xo + Ai + A2 + C2(AoAi + AqAs + A1A2) + C3A0A1A2), 



with 

C2 = 1 + C2, C3 = C2 + C3, (B.9) 

and 

Ao = 1 - e-"/^ (B.IO) 

Ai = 1 - _ ;l - e-^/2(l - A2 - rA'2), (B.ll) 

A2 = l-e-^/2^ (B.12) 

Using expressions (B.7,B.8,B.9) and expression (B.IO, B.ll) in eqs.(B.4,B.5,B.6) gives the 
following system of three equations for three unknowns A, and A2: 

1 - e-^ + re-^X' = -mV2(2A2 + C2XI + (1 + 2C2A2 + C3A^)(1 - e-^/\l - A2 - rA'2))) (B.13) 

1 - e-^ - re-V = -mV'(2A2 + csA^ + (1 + 2C2A2 + C3A^)(1 - e'"/^)) (B.14) 

riy'{l + 2C2A2 + C3XI) = 4(e^/2 - 1)(1 + C2A2 + (cs + C3A2)(1 - e""/")). (B.15) 
After hnearizing in u, A, equations (B.13, B.14, B.15) can be solved exactly: 




(B.16) 



manifesting the vDVZ discontinuity. When mrg <S 1, we have c = 4rg/3 where rg — M/Ai:. 

Without hnearizing, equations (B.13, B.14, B.15) cannot be solved as written. Although 
we do have three equations for the three unknowns, it is seen that the only derivatives of the 
unknowns in these equations are A' and Ag in eq.(B.13), v' in eq.(B.14), and v' in eq.(B.15). 
What makes it possible to solve the system is to derive an algebraic relation between r, 
I/, A, A2 by equating the two expressions for v' . Then one can select any two of the three 
unknowns A, A2 and derive a system of two first-order differential equations for the two 
selected unknowns. 

However, even in the current form one can use Vainshtein approximation to extract a 
solution which is valid well inside the Vainshtein radius r <^ ry = (r^/m^)^/^, but before 
getting too close to the horizon, where Aq diverges as 

(1 -rg/r)-V2 and Tj^ ~ m Ao can no 

longer be neglected from the Einstein equations. 

When ^ r ^ ry, we substitute the linearized Schwarzschild solution 1/ = —X~ —fg/r 
into (B.15), also linearized in v and A, but exact in A2, and solve for A2. We find 

A^ = l/c3. (B.17) 
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which imphes C3 > 0. Numerical integration shows that the positive root is connected to the 
asymptotically decaying solution. 

When r — )■ but still Gm'^{l — rg/r)~^^'^ <^ 1, and assuming that Ai and A2 remain finite 
we expect the geometry to be close to Schwarzschild. Substituting e'^ — e"^ = 1 — Vg/r in 
(B.15) and solving for A2(r) gives, at r = r^,^^: 

A2 = -2 - C2/C3 + ^4 + cl/cl - I/C3 , (B.18) 
r,(l - r,/r)V2A'2 = -A2 - 2C2/ V4ci + ci - C3 . (B.19) 

The same substitution e"^ = 1 — Vg/r in (B.ll) gives, at r = r^, 

Ai = l + rg(l-r,/r)i/2A'2. (B.20) 

These a posteriori justify our assumption of Ai,2 remaining finite. In the limit m — )> 
these approximate solutions are expected to become exact, which is confirmed by numerical 
integration. 



B.l Stress-energy components 

The energy density and stresses in various regions of interest can be obtained from the above 
approximate solutions: 

When <^ r <S m'^, it is seen from eqs. (B.ll), (B.12), and (B.16) that A2 jJi/'i and 
Ai ~ A2 + rAa are the largest eigenvalues, so we can approximate 

r° = m20(r,/r), (B.21) 
r/~-2T|~-^m2f^)'. (B.22) 



3 \r J 

When Tg <^r <^rv, still Ai ~ A2 I/a/^s are the larges eigenvalues, and we have 

T° ~ -m\A + 3c2/V^)/v/^, (B.23) 
Tl ~ T| ~ -m\2 + C2/V^)/V^. (B.24) 

Finally, when r ^ Vg and Aq — > — 00 

T° = m^O{l), (B.25) 
Tl ~ -m^Aoll + 2C2A2 + C3A2), (B.26) 
T| ~ -m'Ao(l + C2(Ai + A2) + C3A1A2), (B.27) 

where Ai^2 are given in (B.18) and (B.20). Pressure p = — (T/ + 2T|)/3 diverges at horizon 
unless C2 = C3. Note that since the energy density Tq is negative on average, p —00 at 
horizon should not be very surprising - reversing the sign of the action does not change the 
equations of motion, but it reverses the sign of stress-energy tensor. 



^^The larger root was chosen for A2 because in the "Vainshtein region" A2 = 1/ ^/c3 which is larger than 
both roots in the parameter range of interest < C3. 
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C Near- hole solutions in FP2 



As argued in §3.3.2 the quasi-stationary near-hole solutions can be obtained from V^T^ = 0, 
or equivalently, from the equations of motion for 0" on Schwarzschild geometry. To derive 
these equations, we first need to determine the eigenvalues \a of the matrix 5^ — ^/H^-, in 
terms of which the FP action is defined. For the spherically symmetric metric 



ds'^ = CdT"^ - 2DdTdR - AdR^ - BR'^dQ^ 



they are given by 



A2 



1 - 
A3 



1 



/2A 

1 - 



A + C± ^/{A^^Cy^^AD^j , 



where 



(C.l) 

(C.2) 
(C.3) 

(C.4) 



A = AC + 

and we have renamed Aq = A+ and Ai = A_. Since the FP2 action (23) depends on A± 
through the combinations A+ -|- A_ and A_|_A_, it is useful to define 



V2 

in terms of which 



1/2 



A + C+^/{A-Cy-4:D^ + A + C- ^/{A - Cy - 4:D^ 



1/2 



A+ + A_ 
A^A- 



= 2-F/VA, 

= i + (i-f)/Va. 



The FP2 Lagrangian becomes 



F 

U ^ 



A 

where we defined the new parameters 

C2 = 1 C2, 

and the functions 

13 = 1 + 2C2A2 + csA^, 



;5 + /? + 2A2 + CaA^, 



C3 = 62 63, 



13 = 1 + 202X2 + csA^ 



(C.5) 



(C.6) 
(C.7) 



(C.8) 



(C.9) 



(C.IO) 



The stress-energy tensor components can now be calculated by varying the action with 
respect to the metric: 



■"^^One can see from these expressions why oc SSpp/SD (and similarly Tqi) is proportional to D with a 
regular proportionality coefficient at D = 0. Spp depends on D via X±, and in the overall ^—g. Both 
are quadratic in D except when: (i) A = C, in which case since A+ = A_ at D = the symmetry of Sfp 
in A± ensures that terms linear in D cancel, (ii) AC = which is a singular point at = and therefore 
cannot be approached from D = side, which is dictated by the asymptotic condition. 
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For instance, we get 



m To 



C 



dU_ 
dD 



dA VAF 



(C.12) 



After deriving the stress-energy components one can use the Vainshtein approximation 
(53), namely. 



C=l-rg/r, A = C-\'y-^-D^), 



(C.13) 



to express everything in terms of r and the two unknowns A2 and D. In particular, using 
R — B~^/'^r, we get 

^ = 7 = di?/dr = 1 - A2 - rA'2. (C.14) 



(C.15) 



The equations for D and A2 can be derived from the stress-energy conservation 



which is particularly simple for the z/ = component in the quasi-stationary approximation. 
It becomes the condition of steady energy flux, 9i(^/^Tq) = 0, whose integration constant 



is (minus) the accretion rate fg = — am r^. Using the expression (C.12), we obtain 



DP 
F 



BR^ 



ar„ 



(C.16) 



Another independent equation can be obtained from the u = 1 component of (C.15) which 
together with the z/ = component form a linear system of equations for D' and X'.^. However, 
it is technically easier to derive this second equation by using the Vainshtein approximation 
directly inside the FP2 Lagrangian (C.8) and writing it as an action for the scalar fields 0" 
on Schwarzschild background: 

S= J drr'^U. (C.17) 

Here U is a function only of r, D, A2, and 7. As a consistency check, note that since 
0°' = D'y / C the equation of motion for 0° is given by 



5S_ 
5^ 



— dr 



7 



— dr 



0, 



(C.18) 



which integrates to our first equation (C.16). The equation of motion for 0* = (1 — A2)rn* 
can be obtained by varying S with respect to A2, but remembering that (fP also depends on 
A2 via 7 = 1 — (rA2)'. Subtracting that contribution, we get 



0X0 



dU DdU 



97 7 dD 



= 0, 



which after some algebra, and using the identity 

CF2 + L>2 = 7-2(1 + 7C)^ 



(C.19) 



(C.20) 
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yields 



rdr 



/3r^(l + 7C) 



(C.21) 



This equation contains both and D' , but it becomes algebraic in A2 once we eliminate D 
using 



2^4 X -1/2 



(C.22) 



which can be obtained by substituting D = ar'^F/r'^P from (C.16) into the identity (C.20). 
After this substitution, the left hand side of (C.21) which would have contained A2 and D' 
becomes 



rdr 



(C.23) 



It is manifestly independent of A2. Factors of A2 will also cancel from the two sides, and we 
finally get 



- 5) + (C2 + C3A2)(1 + (1 - A2)C) = ( C + 



2 4\ 1/2 
2r-i 



1 + C2A2 + 2(C2 + C3A2)). (C.24) 



For given a and r^, A2 is solved for from (C.24), using which D is given by 

.2^4 X -1/2 



K7-' + c)(c + 2^) 

which is obtained by substituting (C.22) in (C.16). 



(C.25) 



D Transition region 

It is difficult to derive analytic results for the transition region r ~ ry = (r^ /m2)i/3, but 
since accretion effects rapidly disappear away from the horizon one can safely linearize in D 
and use the knowledge about static solution in this region. After this linearization, the only 
component of the stress tensor that depends on D is (c.f. eq.(C.12)) 

~ m'^-^S-—, (D.l) 
° %/A(v^+v^)' ^ ^ 



the rest being given by their a = expressions. Here A ~ AC and the Vainshtein approxi- 
mation (C.13) is not used, so (D.l) is valid even for r > ry. 

Inside the Vainshtein radius r <^ ry, we have, from eq.(C.16), a steady flow with energy 
flux 

^2 



ATo^ = am^-f, (D.2) 
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corresponding to Vg = —am'^rg, which in turn fixes the hnearized solution (43) of §3.3.1 and 
from (D.l) we get for ry <^ r ^ m^^ 

d 2 t"^ 

y/Kr^ ~ m^- = -am^^. (D.3) 
2 3 v 

So there is a net positive fiux of energy into the transition region. This can be understood 
in the following way. On a static solution, the relevant degree of freedom in this region is 
R/r = or the related quantity Aa (see eqs.(B.12), (B.16), (B.17)) 

which is approximately a function only of the ratio r/ry- The energy density profile is 
therefore a function of r/ry] as was shown in appendix B.l, in going from r ^ ry to r <^ ry 
it interpolates between a negligible amount to a negative constant e ~ — m^. Once we 
allow accretion, gradually decreases and as a result the Vainshtein radius and the above- 
mentioned profile adiabatically shrink to smaller radii. Therefore, the total amount of energy 
in this region AM ~ evy ~ —Vg grows, requiring a net positive energy fiux. 

The exact value of AM can be found from the difference of the invariant masses Mi 
and M2 inside two spheres of radii ri and r2, where ^ ri ^ ry and ry <^ r2 ^ m~^. 
Since at distances r ^ Vg the geometry is nearly fiat, the metric can be transformed into 
dfiu = V^iv + h^i, with h^y ^ 1, in terms of which the invariant mass inside a sphere of radius 
r is given by 

M = ^ y {d.hjj - djhij)n'r^dn. (D.5) 

Using circumference defined radius to parametrize the static metric [as in (B.2)], we get hij = 
—Xn^n\ where inside the Vainshtein radius A is given by the Einsteinian Vg/r expression, 
and outside by the vDVZ one 2rg/3r. The enclosed masses are then 

Ml = inrg, M2 = yr^, (D.6) 

giving 

AM = M2-Mi = -—rg. (D.7) 

Therefore, when a ^ the time-dependence of AM exactly accounts for the difference 
between (D.2) and (D.3). 

Similar methods can be used to find D throughout the transition region. One uses 
the numerical static solution of [17] to determine M{r,rg). By numerically differentiating 
M{r,rg) with respect to Vg, one can find 

M{r,rg) = rgdr^M{r,rg). (D.8) 

The energy fiux at any r2 ~ ry is determined from 

M{r2) - Ml = -r^^Tl , (D.9) 

using which one solves for D from (D.l). 
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E The mass of a star and the global charge Qq 



Consider a star of gravitational radius Vg deep inside the Vainshtein regime. In the Schwarzschild 
coordinates, where the metric is given by 

ds'' = e'^dt' - e^dr'' - r^dQ'', (E.l) 

we have e" = = 1 — Vg/r for r <^ ry- We will show in what follows that for finite redshift 
stars the Noether charge Qq coincides with minus the Schwarzschild mass Ai^Vg of the star. 

Note first that for all r ^ Vg the functions v and A are infinitesimal (although different 
from the Einsteinian value when r » ry), so that the space-time is nearly flat. Thus, the 
ADM mass within any sphere of radius r ^ is well-defined. Consider a sphere of radius 
Ti where Tg <^ ri <^ ry. The enclosed ADM mass coincides with the Schwarzschild mass 
-i-RTg. On the other hand the ADM mass inside a sphere of much larger radius, r ^ m~^, 
must vanish in massive gravity; therefore, the contribution of Spp to Madm coming from 
the region ri < r < oo should be exactly —AnVg. Given the closeness of the metric (E.l) 
to Minkowski in this region, this contribution of Spp can be directly calculated from the 
volume integral of its energy density, namely 



d\Tl = -Anrg. (E.2) 



We will show that the integral on the left coincides with the volume integral of Jq in the 
same region of space, which is what one expects if Qo is a good characterization of the energy 
of Spp- 

Since for r ^ Vg the unitary-gauge metric is approximately diagonal (even if there is 
time-dependence), we have (p^ c:^ t and (9^0* ^ 0, so Jq and can be replaced with the 
unitary-gauge expressions Jq and Tq . From (76) we have (after redefining m'^U U) 

Jo = ^0 - U, (E.3) 
and the above requirement on Qo implies 

drr^U < r-g. (E.4) 



/ 



'riS>rg 

This is trivially satisfied in the linear regime r ^ ry, where the perturbations of the unitary- 
gauge metric hab are themselves infinitesimal, and U — 0{m'^hl^). 

Now consider the contribution from a region ri < r < r2 where <C ri <C ry <S r2 <^ 
m~^. The static solution is a good approximation throughout this region and we can use 
the formalism of appendix B to verify eq.(E.4). Since u <^1, the stress-energy conservation 
(B.6) simplifies to 

Tl'c,^{Ti-Tl), (E.5) 

and for a diagonal unitary-gauge metric, and T2 of a generic FP (including FP2) are 
given by 

T; = (X.-1)§-^-U, r| = l(A.-l)||-f/. (E.6) 
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(the factor of 1/2 in T2 is because A2 = A3). Moreover, in this region A ^ A2, and the 
expression (B.ll) for Ai can be approximated as 



Ai ^ (Aar)', 



(E.7) 



while Aq 
yields 



Z//2 and its derivatives can be neglected. Substituting (E.6) and (E.7) in (E.5) 



We can now transform 

fr-2 



drr U = -r^U 



ri 



r2 



dU 
1 

~ 3 
1 

" 3 



(E.8) 



r2 



drr^U' 



ri 
r2 



drr 



ri 



dX 



9Ai 



and use (E.8) to obtain 



r2 



drr'^U 



U - rX' 



dU 
dX[ 



r2 



ri 



(E.9) 



(E.IO) 



But r2 is in the linear regime where Xa ~ haa <^ 1, and hence the right hand side which is 
quadratic in A^ is negligible at the upper limit. At the lower limit, since U oc m^, Ai 2 ~ 1, 
and rf Ty = Vg/m'^, we get the desired relation (E.4). 

For a star of finite redshift all eigenvalues A^ remain finite and the contribution of the 
region r < ri to Qq, which is of order rjjg ~ m^rf, is much less than Vg. Therefore, Qq is 
given by the total amount of negative mass in the region r ^ Vg evaluated by the volume 
integral of Tq and we have Qo = —Airrg. 



F Fluid accretion 

Potential flow of ideal fluids can be described by a single scalar fleld (p with the k-essence 
action [34] 

S = J d'x^gP{X), (F.l) 

where X = g^^dfj,(f)d,j(j). The stress-energy tensor is 

V = '^Pxd^<pd,<P - g^^P, (F.2) 

where Px = dP/dX. We can therefore identify the fluid pressure p, energy density e, and 
velocity as 

p = P, e = 2PxX-P, u^ = d^<j)/Vx. (F.3) 
The p = c^e equation of state, therefore, corresponds to 

1 + 

P = X", with n = ' . (F.4) 
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3.5 4.0 



Figure 2: f{r,(p') plotted as a function of if' at fixed values of r and for rg = 1. When < 1 
(left), the maximum of /(r, (p') decreases until a critical radius is reached and then increases. When 
> 1 (right), f{r,ip') has no maximum. 



For static spherically symmetric fluids with asymptotic density sq, we have (j) = (pot, with 
00 = {c^SoY^^^'^^'' \ which on Schwarzschild geometry gives the singular at horizon solution 
(61). To have an accreting solution we take (p = (j)o{t + (p{r)). Integrating once the (p equation 
of motion, we get 

/(r, ^ r\l - "fMil - ^)-i - (1 - ^)y,'2]n-i = a, (F.5) 

where A is an integration constant. For any A, there exists at large values of r a decaying 
solution if' ^ A/r'^, corresponding to the steady accretion 

T; = -{l + cl)eo^^. (F.6) 

However, when < 1 only one of the decaying solutions matches a regular solution at 
horizon (the Bondi solution). This can be shown as follows. At any fixed r the function 
f{r,(f') in (F.5) is maximized at 

^' = c.(l - ^)-\ (F.7) 
r 

giving a maximum fma.x{r)- As r is decreased from infinity, fma.x{r) reaches a minimum at 
the critical radius 

rc = -^^rg, (F.8) 

and then increases (figure 2). For A > ^Bondi = /maxl'^c), the asymptotically decaying 
solution ceases to exist below some radius r(> r^). For A < A^ondi, there exist a solution all 
the way to the horizon but (p' is not large enough there. Therefore, X diverges as (1— r^/r)"^, 
leading to a pressure singularity. Only for A = ^Bondi the asymptotically decaying solution 
connects to a rapidly falling solution at horizon, where v?' — ?■ (1 — Tg/r)^^ and X remains 
finite. 
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When > 1, f{r,ip') is unbounded above and there is no critical point as can be seen 
from figure 2. For all A the asymptotically decaying solution goes into 



^' = (1 - -^r - o > (F-9) 



2 

in the r ^ Vg limit, and has a finite pressure. The energy of fiuid elements at horizon reaches 
the constant 

References 

M. S. Volkov, Phys. Rev. D 85, 124043 (2012) [arXiv:1202.6682 [hep-th]]. 
R. A. Konoplya and C. Mohna, Phys. Rev. D 75, 084004 (2007) [gr-qc/0602047]. 
M. Fierz and W. Pauli, Proc. Roy. Soc. Lond. A 173, 211 (1939). 
D. G. Boulware and S. Deser, Phys. Rev. D 6, 3368 (1972). 

C. de Rham and G. Gabadadze, Phys. Rev. D 82, 044020 (2010) [arXiv: 1007.0443 
[hep-th]]. 

C. de Rham, G. Gabadadze, A. J. ToUey, Phys. Rev. Lett. 106, 231101 (2011). 
[arXiv:1011.1232 [hep-th]]. 

T. M. Nieuwenhuizen, Phys. Rev. D 84, 024038 (2011) [arXiv:1103.5912 [gr-qc]]. 

S. F. Hassan and R. A. Rosen, Phys. Rev. Lett. 108, 041101 (2012) [arXiv: 1106.3344 
[hcp-th]]. 

S. F. Hassan and R. A. Rosen, JHEP 1204, 123 (2012) [arXiv:1111.2070 [hep-th]]. 

M. Mirbabayi, Phys. Rev. D 86, 084006 (2012) [arXiv:1112.1435 [hep-th]]. 

C. de Rham, G. Gabadadze, L. Heisenberg, D. Pirtskhalava, Phys. Rev. D83, 103516 
(2011). [arXiv:1010.1780 [hep-th]]. 

A. Gruzinov, [arXiv: 1106. 3972 [hep-th]]. 

C. de Rham, G. Gabadadze, A. J. ToUey, [arXiv: 1107.0710 [hep-th]]. 

S. Deser and A. Waldron, Phys. Rev. Lett. 110, 111101 (2013) [arXiv: 1212.5835 [hep- 
th]]. 

L. Berezhiani, G. Chkareuli, and G. Gabadadze, unpublished. 
N. Arkani-Hamed, H. Georgi and M. D. Schwartz, Annals Phys. 305, 96 (2003). 



31 



[15] S. L. Dubovsky, JHEP 0410, 076 (2004). [hep-th/0409124]. 

[16] A. H. Chamseddine and V. Mukhanov, JHEP 1008, Oil (2010) [arXiv:1002.3877 [hep- 
th]]. 

[17] A. Gruzinov and M. Mirbabayi, Phys. Rev. D 84, 124019 (2011) [arXiv:1106.2551 [hep- 
th]]. 

[18] A. Salam and J. A. Strathdee, Phys. Rev. D 16, 2668 (1977). 

[19] K. Koyama, G. Niz, G. Tasinato, Phys. Rev. Lett. 107, 131101 (2011). [arXiv: 1103.4708 
[hep-th]], 

K. Koyama, G. Niz, G. Tasinato, Phys. Rev. D84, 064033 (2011). [arXiv:1104.2143 
[hep-th]]. 

[20] L. Berezhiani, G. Chkareuh, C. de Rham, G. Gabadadze and A. J. ToUey, Phys. Rev. 
D 85, 044024 (2012) [arXiv:1111.3613 [hep-th]]. 

[21] C. Deffayet and T. Jacobson, Class. Quant. Grav. 29, 065009 (2012) [arXiv: 1107.4978 
[gr-qc]]. 

[22] T. Jacobson, Phys. Rev. Lett. 83, 2699 (1999) [astro-ph/9905303]. 
[23] See page 21 of 

http:/ / www.yukawa.kyoto-u.ac.jp/ contents / seminar / archive/2003 / str2003 / talks / dvali.pdf 

[24] H. van Dam and M. J. G. Veltman, Nucl. Phys. B 22, 397 (1970); 

V. L Zakharov, JETP Lett. 12, 312 (1970) [Pisma Zh. Eksp. Teor. Fiz. 12, 447 (1970)]. 

[25] L. Hui and A. Nicolis, arXiv: 1202. 1296 [hep-th]. 

[26] L. Hui and A. Nicolis, Phys. Rev. Lett. 109, 051304 (2012) [arXiv: 1201. 1508 [astro- 
ph.CO]]. 

[27] A. L Vainshtein, Phys. Lett. B 39, 393 (1972). 

[28] G. Chkareuh and D. Pirtskhalava, Phys. Lett. B 713, 99 (2012) [arXiv: 1105. 1783 [hep- 
th]]. 

[29] E. Babichev, C. Deffayet and R. Ziour, Phys. Rev. Lett. 103, 201102 (2009) 
[arXiv:0907.4103 [gr-qc]]. 

E. Babichev, C. Deffayet, R. Ziour, Phys. Rev. D82, 104008 (2010). [arXiv: 1007.4506 
[gr-qc]]. 

[30] L. Berezhiani, G. Chkareuh and G. Gabadadze, arXiv: 1302.0549 [hep-th]. 
[31] H. Bondi, Mon. Not. Roy. Astron. Soc. 112, 195 (1952). 

[32] E. Babichev, V. Dokuchaev and Y. .Eroshenko, Phys. Rev. Lett. 93, 021102 (2004) 
[gr-qc/0402089]. 

32 



A. D. Dolgov, H. Maeda and T. Torii, hep-ph/0210267. 

C. Armendariz-Picon, T. Damour and V. F. Mukhanov, Phys. Lett. B 458, 209 (1999) 
[hep-th/9904075]. 



33 



